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The collinear color electric and magnetic fields have been discussed to be produced immediately 
after high energy heavy ion collisions. We discuss anomalous gluon production under the background 
gauge fields. The gluons are Nielsen- Olesen unstable modes. The production rate of the modes by 
Schwinger mechanism has recently been found to be anomalously larger than the rate of quarks or 
other stable gluons. Analyzing classical evolutions of the modes with initial conditions given by 
vacuum fluctuations, we find that their production makes the color electric field decay very rapidly. 
The life time of the field is approximately given by the inverse of saturation momentum in the 
collisions. We discuss a possibility that the mode with zero momentum form a Bose condensate and 
its gluon number density grows to be of the order of l/a g - The condensate would melts into quark 
gluon plasma after the gluon number density grows sufficiently large for nonlinear interactions to 
be operative. 
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I. INTRODUCTION 



It has recently been paid attention to that color electric ( E ) and magnetic ( B ) fields produced in high energy 
heavy ion collisions [lj decays by producing quarks and gluons. Namely, glasma decays into quark gluon plasma ( 
QGP ). One of the most efficient driving forces for the decay is the Schwinger mechanism The production rate 
, of massless quarks by the mechanism is similar to the production rate of gluons. The difference simply comes from 
the difference in the color charges, spins and statistics. However, when the color magnetic field ( B ) is present, the 
\ production rate of the gluons becomes anomalously larger than the production rate of the quarks. This important 
qq ■ fact has recently been shown by Itakura and TanjijS]. They have found that the production rate of the gluons is 

proportional to exp(i:B/E), while the rate of the quarks is proportional to exp(—TtB/E) or unity. 
C\| . As is well known, the gluons can have imaginary energies under the color magnetic field. Namely, the states 
with magnetic moments parallel to B in the lowest Landau level are unstable; their energies ( oc — gB ) are 
imaginary when p| < gB, where P3 is a momentum component in the direction parallel to B — (0,0, B) ( g is 
a gauge coupling constant ). Thus, their amplitudes exponentially grow with time. The amplitudes of the other 
states of the gluons simply oscillate with time, because their energies are real. The presence of such unstable modes 
( Nielsen-Olesen unstable modes[6J ) has prevented us to find a proper quantization of the gluons. Hence, the 
production of the Nielsen-Olesen ( NO ) unstable gluons by the Schwinger mechanism was not discussed. But the 
above authors pointed out that the unstable modes become stable owing to the acceleration by the electric field. This 
is because the momentum ps(t) = f dt' gE becomes large sufficiently such that the imaginary energy becomes real 
with Ps(t) > (gB) 2 . Consequently they have derived the production rate of the NO gluons. 

The importance of their result is that the production rate, in other words, the occupation number ( oc exjp(nB / E) 
) of a gluon state associated with the NO unstable gluons is much larger than the unity. This should be compared 
with the typical occupation number ( oc exp(—irB/E) ) associated with other stable gluons or massless quarks [3]- It 
is much less than the unity. The large production rate of the gluons leads to the rapid decay of the electric field. In 
particular, the decay process accelerates because the production rate oc exp(-7r B/E) increases more as the electric field 
E becomes weaker. Hence, when we discuss the decay, it is important to include the back reaction of the electric field 
to the production of the gluons. Furthermore, the fact of the large occupation number allows us to classically analyze 
the production of the NO gluons. When we classically analyze the production, we can easily take into account the 
back reaction of the electric field to the gluon production. The back reaction has not yet been discussed by the above 
authors 0. 

In this paper, by assuming that the NO modes are initially produced by vacuum fluctuations, we analyze the 
classical evolution of the modes and the color electric field. We find that the electric field rapidly decays owing to 
the acceleration of the large amount of the NO gluons. The life time of the electric field is approximately given by 
QjT 1 where Q s denotes saturation momentum in high energy heavy ion collisions. We show that the color electric 
current carried by the gluons is much larger than that of the quarks produced by the Scwhinger mechanism. Thus, the 
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influence of the quarks on the decay of the electric field is negligible compared with that of the gluons. We also discuss 
that the NO gluons form a Bose condensate, because of the large occupation number of the gluons. Especially, the 
amplitude of the NO mode with zero momentum is the largest and grows most rapidly. It can grow to be of the order 
of 1/g until nonlinear interactions of gluons are operative. Then, it would melt into QGP through the equipartition 
of the momentum by the nonlinear interactions. It implies that the number density of the gluon condensate grows up 
to of the order of 1/g 2 - 

In the next section [TT] we briefly review the Nielsen- Olesen unstable modes and discuss anomalous production of 
the modes. In the section IIII1 we discuss field configurations of the Nielsen- Olesen unstable modes and find basic 
equations governing the temporal behaviors of the modes as well as the color electric field. In the section IIV1 we find 
that the gluons of the modes are dominantly produced and their production leads to the rapid decays of the electric 
field. In the section[V]we discuss the gluon condensation of the NO modes, which may arise before the thermalization 
of the quark gluon plasma. In the final section IVil we summarize our results. 

II. NIELSEN-OLESEN UNSTABLE MODES 

We first explain our formalism and briefly review the NO unstable modes. We also explain the anomalous production 
of the NO modes under the electric field. We consider SU(2) gauge theory and suppose the background color electric 
and magnetic fields given by E a = S a ^(0, 0, E) and B a = 6 a ,3(0, 0, B). They are supposed to be spatially homogeneous 
and collinear both in the real and color spaces. The gauge fields are represented by the gauge potential = A^ 3 . 
Under the background fields, the gauge potentials $ M = (A 1 +iA 2 ^) / 'v2 perpendicular to A 3 ^ behave as charged vector 
fields. When we represent SU(2) gauge potentials A^ using the variables A^ and (f^, Lagrangian of SU(2) gauge 
potentials is written in the following, 

L = -\fI„ \\D^ v - D v *tf - ig(d»A v - d v A M )&»$" + £@t$u <^) 2 , (1) 

with Ffj, >u — d^Av — d u Afj, and = — igA^, where we used a gauge = 0. Obviously, we find that the 

fields $ M represent charged vector fields with the anomalous magnetic moment described by the term —ig{d^A v — 
dvA^Q^Q" . Therefore, it is easy to see that when the magnetic field B = d\Ai — &2A\ is present, but E = 0, 
the particles represented by the fields $ occupy the Landau levels and interact with each other through the term 
^j-("I > j J < I>„ — $J,<E> A1 ) 2 . The energies of the states in the Landau levels denoted by integer N > are given by En = 
y / 2gB(N + 1/2) ± 2gB + p|, where ± denotes magnetic moment parallel ( — ) or anti-parallel ( + ) to B. 

Among them we notice the states in the lowest Landau level ( N = ) with the magnetic moment parallel to B. Their 
energies can be imaginary; -Ejv=o = \Zp§ gB. Thus, the modes with the imaginary energies exponentially increase 
or decrease with time. That is, the field $ representing the modes evolves with time such that $ oc exp(— iE^^t) = 
exp(±|£ , AT =0 |i) = exp(±| \J gB — pf| t). The states are called as Nielsen- Olesen unstable modes. Among them, the 
mode with P3 = increases or decreases most rapidly. The presence of such unstable modes implies the instability of 
the vacuum state, that is, the state with the background color magnetic field B. This is similar to the case that the 
state ijj = is unstable in a model of a complex scalar field with the double well potential — m 2 |-0| 2 + ■||'0| 4 ( m 2 > 

). In this model unstable modes around the state ■0 = exist and exponentially grow such as ip oc exp(ty/ m 2 — p 2 ) 
where p denotes a momentum. Similarly the background gauge fields with color magentic fields are unstable. Indeed, 
the classical simulations have been performed to show the instability of the background gauge fields involving color 
electric and magnetic fields. The instability in the simulation has been discussed to be caused by the NO unstable 
modes: Their amplitudes exponentially grow and then saturate when nonlinear interactions are operative owing to 
the growth of the amplitudes. 

The spontaneous production[l(| of the unstable modes ( or gluons ) is caused by the color magnetic field. Thus, it 
is not the Schwinger mechanism. When we analyze the Schwinger mechanism of the NO modes, the presence of the 
imaginary energy was an obstacle because we cannot properly quantize the modes. But, we should note that in order 
to obtain the production rate by the Schwinger mechanism [3|, we only need in-state ( in infinite past ) and out-state 
( in infinite future ) under the electric field. Although the NO modes are unstable, they are stable in the infinite 
past and future in the presence of the electric field. This is because the square of the momentum p 3 (t) — J dt' gE is 
sufficiently large in the past and future such as p 2 (t) — gB > 0. The stability of the NO modes in the infinite past 
and future allows us to estimate the production rate[5| of the modes under the electric field. 

In general, charged fields oscillate with the frequency cx gEt in the past ( t — > — oo ) and future ( t — > +oo ) under 
the electric field E, owing to the acceleration by the electric field. The frequency is real and depends on time. In 
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this case the production rate of the fields is less than unity; it is proportional to, for example, exp(—TrB/E) when 
the magnetic field is present. But, in the case of the NO modes they pass a period in which their frequency becomes 
imaginary so that their amplitudes exponentially grow. Before and after the period they simply oscillate with the real 
frequency. In other words, the modes oscillate and their amplitudes smoothly change in the far past, but once they 
enter the unstable period, their amplitudes exponentially grow with time. After passing the period, they oscillate 
again but with much larger amplitudes than those before passing the period. These behaviors are peculiar to the NO 
modes. In particular, the exponential growQ of the amplitudes leads to the anomalous production rate oc exp(nB/E). 

We may understand a naive physical reason why the rate increases more as the electric field becomes weaker. The 
large production rate comes from the fact that the modes pass the period in which they exponentially grow. They 
stays the period approximately for the time At = \J gB j gE because the momentum increases within the time such 
that Ap 3 = AtgE = y 1 gB. Hence, the amplitude grows by exp(^gBAt) = exp(B/E). Weaker electric field causes 
the longer stay in the unstable period and hence the larger grow of the amplitude. Although it is a rough estimation, 
it explains why the production rate of the gluons oc exp(7r B / E) becomes large as the electric field becomes weak. 

Hereafter we take only the NO modes and analyze their production under the electric field. The NO modes are 
described by the field <I> = ($1 + i<i>2)/v / 2 and are governed by the following Hamiltonian 

H = J Sx{^{d A 3 fMd^\ 2 +\(d^gA)^-2gB\^) = J d 3 x(^(d A 3 ) 2 + \d <S>\ 2 + \(d 3 -tgA 3 )<S>\ 2 -gB\<i>\ 2 ) , (2) 

where we neglected the nonlinear interactions ^($>j^$ u — ^"I^) 2 in eq([T]). The color electric and magnetic fields are 

given such that E = OqA = (0,0, 3qA 3 ) and B = d x A — (0,0, B). The nonlinear interactions are not operative as 
long as the amplitude <!> is small. When the field $ grow large such as the nonlinear interactions are operative, the 
NO modes couple with the other modes in higher Landau levels as well as themselves. 

We note that the last term in the Hamiltonian represents a negative potential. When the magnetic field forms a 
flux tube, the term represents a negative potential with finite width given by the width of the flux tube. Whether or 
not the field <I> possesses unstable modes depends on the existence of the states trapped in the negative potential. 

It is interesting to see the analogy between our model and the model of the complex scalar field with the double 
well potential mentioned above. The model describes Cooper pairs condensates. Thus, the decay of the color electric 
field corresponds to the decay of an external electric field imposed on superconductors. The question is how fast the 
external electric field decays, just after a metal in a normal state is supercooled with a temperature below the critical 
temperature separating superconducting ( (ip) = \J m 2 /\ ) and the normal states ( ("0) = ) . The normal state ip — 
decays by producing the Cooper pairs ip, which condense to form the state (ip) — y/ m 2 /A. Since they are accelerated 
by the electric field, the electric field loses its energy and vanishes. 

III. PRODUCTION OF NIELSEN-OLESEN UNSTABLE MODES 

Because the production rate is much larger than unity ( this implies that the occupation number in a state is much 
larger than unity ), the production of the NO gluons represented by NO modes may be classically analyzed. Then, 
we can easily take into account the back reaction of the electric field to the gluon production. In this section we will 
formulate basic equations governing the back reaction. 

First, we discuss the assumption that the background gauge fields are homogeneous in the transverse plane per- 
pendicular to the collinear fields B and E. When the NO unstable modes are excited, they destroy the homogeneity 
because of the localization of the wave functions of the modes. 

gB\z\ 2 

= (xx - ix 2 ) n exp( - h ip 3 x 3 ), (3) 

with z = x\ +1x2 and integer n > where we used a gauge potential A = (—BX2/2, Bxi/2, 0). The effect of the back 
reaction induces the inhomogeneity in the electric field; the currents carried by the modes are not homogeneous so 
that the electric field affected by the currents is also not homogeneous. 

But, by taking the appropriate linear combination of the NO unstable modes we can form almost homogeneous field 
configurations in the transverse plane. Then, their currents are also almost homogeneous. Such field configurations 
are given by, 
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* = H <M f )> "M^ = / dp 3 c{p 3 )exp( +ip 3 x 3 ), (4) 

with c{p 3 ) is a dimensionless function of the longitudinal momentum p% and Z{ — x\.i + 1x2,%, where each component 
(j>i satisfies the condition, (j>i4>j — $i j4>i since we impose that \zi — zA > Ib = , ■ Namely, a configuration ^ is 

y/gB 

separated with the nearest neighbors approximately by the distance Ib- Furthermore, we assume that the area L 2 of 
the transverse plane is given by L 2 = Nl 2 B . Thus, we find that the field configuration 4> is approximately uniform in 
the transverse space. This kind of the configuration of the NO modes was analyzed[ll| to discuss so called "spaghetti 
vacuum" . 

Using the field configuration, we rewrite the Hamiltonian of the NO unstable modes, 




dz—gA)<&\ 2 ~gB\<&\ 2 j ~ N J d 3 x\^(d A) 2 M^\ 2 M(id 3 -gA^\ 2 -gB\^\ 2 i (5) 

with <f> = j dp 3 c{p 3 ) exp(ij»3X3 — gB\z\ 2 /4), where the color electric field is given by E — 3qA with the homogeneous 
gauge potential A = A3 . The Hamiltonian describes the NO modes under the homogeneous background electric and 
magnetic fields. The first term represents the energy of the electric field and the other terms represent the energy of 
the NO modes. We can see that the last term with the magnetic field gB represents a negative potential for the NO 
modes <f). Thus, it gives rise to the imaginary energy of the field </> oc exp(— iet) with e 2 = — gB < 0. 

If the magnetic field forms a flux tube with a finite width, it gives a negative potential with the finite width. Thus, 
if the field is trapped by the potential, the energy e can be imaginary, but its absolute value is smaller than \/gB; it 
depends on the width of the tube. When the width is infinite, the energy e is given by y '—gB. On the other hand, 
when the width is finite, as it becomes smaller, the absolute value of energy e becomes smaller. The flux tubes of the 
background gauge fields are more realistic field configurations produced in high energy heavy ion collisions than the 
homogeneous ones under consideration. Since the gauge fields are homogeneous in the longitudinal direction, they 
can be viewed as an ensemble of electric and magnetic flux tubes with various widths. Based on the view, we can 
understand^, Q the results of the numerical simulations @- Although the flux tubes are realistic one, it is meaningful 
to analyze the anomalous gluon production in the homogeneous gauge fields in order to see physical essences of the 
production. 

We proceed to analyze the back reaction. For the purpose, we decompose the field 4>{x) into the components of the 
momentum eigenstate, 

j=^r<f> p exp(ip£) (6) 
v 

with 



8 " 2C(P3) exp(-4) (7) 



rp gBL 3 / 2 ^ y gB 
where the transverse momentum px is given such that p\ =p\ +p 2 - Then, it follows that 



H = L 2 (^(d A) 2 + 1 ]T( Ido^l 2 + |(p 3 + gA)<f> p \ 2 - 




(8) 



where we have used the following formulas, 



d 3 xexp(ipf) = (2tt) 3 (5 3 (p) = L 3 <5 3 . . (9) 
The equations of motions of the fields <j) p and the gauge field A are given by respectively 
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d 2 cf> p = gB<t> p - (p 3 + gA) 2 4> v and Ld 2 A = -J- V(p 3 + gA)^ 2 (10) 

where the second equation represents a Maxwell equation d^E = — J with the current J = l£ Ss(P3 + .9^)l^p| 2 - It 
describes how the electric field changes by the effect of the current J. 

We should point out that the largest amplitude of the NO modes is given by the mode with the vanishing transverse 
and longitudinal momentum, that is, the mode with p? = in eq® and + gA — in eq (|10p . The mode grows 
most rapidly compared with the other modes with p 3 + gA ^ 0. Here, the momentum p% + gA(t) = p 3 + g J Q dt'E(t') 
is the momentum gained by the mode with initial momentum p 3 owing to the acceleration by the electric field E. We 
naively expect from the similarity to the model of the complex scalar field that the mode with zero momentum form 
a stable Bose condensate. But as we discuss later, although the mode form a Bose condensate with its amplitude of 
the order of 1/g, the condensate is unstable. 

In order to solve the equations we need to impose initial conditions. The initial condition of the electric field is 
given such that E(t = 0) = doA(t = 0) = E and A(t = 0) = 0. This corresponds to the fact that we switch on the 
electric field E = Eq at t = 0. In other words, we consider the situation that high energy heavy ion collisions occur 
at t = and the coherent color gauge fields are produced at the instance. 

How should we choose initial conditions of the field <f> ? Before the collisions, there are no color electric and magnetic 
fields. The gluons with small x form color glass condensates in nuclei. Just after the collisions the gluons form the 
coherent gauge fields E and B, but there is no classical field 4>. Thus, we may naively choose the initial conditions 
such that do4>(t = 0) = and 4>{t = 0) = 0. But these initial conditions lead to the trivial result; 4>{t) = for any 
time t > 0. Therefore we need to find an appropriate initial condition. As we explained in the previous section, the 
NO modes are spontaneously generated owing to the instability of the state with the homogeneous magnetic field. 
Thus, it is reasonable to take the initial conditions given by the vacuum fluctuations, 



Mt = 0) = y(<^) and d t <t> p (t = 0)= ^J((M P ) 2 ) (11) 

where <f> p denotes the momentum component of the free massless scalar field <j> with no background fields. The state 
| ) represents the vacuum state (<fi) = 0. Then, the vacuum can be represented by the following wave functionals 




with \p\ = \J p\ + p\ ■ Namely, we assume that the initial conditions are given by the vacuum fluctuations in the 
vacuum without E and B. Then, using the wave functionals, we can calculate the expectation values in eq. (|ll[) . As 

the NO modes are given such that 4> p = ^ exp(-f|) with arbitrary dimensionless function a(p 3 ) ^ 8^c(p 3 ), we 
insert this function into the wave functionals and obtain the distributions of a{pz) and doaiji^) in the vacuum, 



#(a(p 3 )) oc exp ( - 

P3 



Hp 3 )\ 2 T( P3 h 
4ttL^B 



and ^ (doa(p3)) oc exp 



E 

P3 



\d a(p 3 )\ 2 U(p 3 lB) 
AttL{ 9 B) 3 / 2 



(13) 



with Ib — 1/y/gB, where the functions T{p 3 Ib) and U(j>31b) are given by 



T(x)= ^ dy^Tx^ exp(-2i/) and U(x) = f°° dy ^LJ^jr . (14) 
Jo Jo \Jy + x z 

Using the wave functionals we calculate the initial conditions for (f> p or 0(^3). Obviously the distribution of the 
real part of 0(^3), that is, 4> p is identical to that of the imaginary part of 0(753). The expectation values of these 
components can be obtained in the following, 



(a(p3 



fg-BL 



J^dzg(z,p 3 ) = T(p 3 l B ) 



and ((d a(p 3 )) 2 ) = 



J^dzz 2 f(z,p 3 ) = n(gBf/ 2 L 
J^dzf(z,p 3 ) U( P3 l B ) 



(15) 
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with 

I \ i z2t (P3 1 b), j p / \ , z 2 U(p 3 l B ) , na , 

g ( z , P3 ) =exp (-— 7 = m -) and /( Z , a> ) = aq) (- i -^_ z ), (16) 

We rewrite the equations of motion using the variables of the real or imaginary part of 0(^3) instead of <fi p , 

d 2 b(p 3 )=gBb(p 3 )-(p 3 +gA) 2 b(p 3 ) and Ld 2 A = £(p 3 + gA) b(p 3 ) 2 (17) 

where the real function b(p 3 ) represents the real or imaginary part of a(p 3 ). Therefore, the initial conditions are given 
by 



A(i = 0)=0, d Q A(t = 0)=E Q and b(p 3 , t = 0) = J^f^r, d b(p 3 , t = 0) = J n( f f )3/ \ L . (18) 

Because these equations involve the size L of the system, it apparently seems that our results depend on the size 
L. But as we show below, our results, e.g. life time of the electric field, are independent on the size L. In order to 

see the independence, we further rewrite the equations of motion by using the variable a(p 3 ) = , / T ^ P3lB '> b(p 3 ), 

d 2 a(p 3 )=gBa(p 3 )~(p 3 +gA) 2 a( P3 ) and d 2 g A = - 9 -^- / + °° dp 3 iP3 t^f]^ (19) 

47T J^ce 1 {p 3 lB) 



with the initial conditions a(p 3 ,t = 0) = 1 and doa(p 3 ,t = 0) = y Hj§^~y^; where we used the relation J2 P3 = 

'Tk J-00 dP3- Obviously both the equations of motion and initial conditions are independent on the system size L. 
Furthermore, the color electric current J is given by 

g^B [ + °° ( P3 +gA)a(p 3 ) 2 

J=— ; / d P3 7fT, , \ . ( 20 

4?r J_ qc T{p 3 l B ) 

which is also independent on L. The current vanishes at t = because T(x) = T(—x) and A(t = 0) = 0. Using 
these equations of motion, we can discuss the temporal behaviors of the electric field do A and the NO modes a(p 3 ) 
by taking into account the back reaction of the electric field to the production of the modes. 
For the convenience, we write down the field <f> of the NO modes in terms of the variable a, 

* = ir^ exp(¥f) = 7l^ w/^Tfe^ exp( "S } exp(¥a °' (21) 

where a{p 3 ) is a dimensionless real function. We can see that the vacuum fluctuation (j>(t = 0) of the NO modes is of 
the order of unity compared with the background gauge fields ~ 0(1/ 17). 



IV. NUMERICAL RESULTS 



Now we wish to discuss the production of the NO gluons and quarks. Especially, we would like to discuss the 
ratio between the amount of the gluons and that of the quarks produced by the electric field. In order to discuss the 
amounts of the particles we compare the color electric current of the quarks with that of the gluons. We show that 
the amount of the gluons is about a hundred times larger than that of the quarks. As a result the color electric field 
rapidly decays owing to this anomalous gluon production. 

First we derive a relevant equation describing the production of the quarks by the Schwinger mechanism. The 
equation has been previously explored[l|- I would like to explain it shortly. We assume that the quarks are massless 
and they form a SU(2) doublet. Then, the charges of the quarks are given by g/2 and —g/2. Both of them possess 
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their anti-particles with their charges given by — g/2 and g/2, respectively. Therefore, we have four massless fermions; 
a pair of the quarks in a SU(2) doublet and their anti-quarks. They are two fermions with the charge g/2 and two 
fermions with the charge —g/2. Their number densities are identical to each other because a pair of a positive and a 
negative charged fermions is created at the same moment under the electric field. 

It has been recently shown 0, [H[ that the evolution of the number density n q of the massless fermions is governed by 
the chiral anomaly when collinear strong electric and magnetic fields are present. In particular the anomaly equation 
is very efficient when the magnetic field is sufficiently strong such that the particles in the production are only allowed 
to occupy the lowest energy states in the lowest Landau level. Then, the equation of the chiral anomaly is given by 



^ = 4a „„, = 2 feZ«, (22) 

where we assumed the homogeneity of the chiral current in the transverse and longitudinal direction; dJ 5 = 0. The 
equality O^Jq = 4dofi q in eq (l2"2l conies from the fact that all of the four fermions have the positive chirality when 
E parallel to B. This is because the positive ( negative ) charged fermions are accelerated to the direction parallel ( 
antiparallel ) to E and their spins point to the direction ( antiparallel ) parallel to B when they occupy the lowest 
Landau level. 

Obviously, the chiral anomaly describes how the number density n q evolves with time under the effect of the electric 
and magnetic fields. The electric field loses its energy owing to the acceleration of the quarks as well as the gluons. 
Hence, we add the contribution of the quarks to the Maxwell equation. Consequently, the equations describing the 
evolution of the number of the quarks and the NO modes as well as the evolution of the electric field are given by 



a E(t)B a /nR f + °° 

d n q = ' o ' d l a (Pz) = gBa(p 3 y(p 3 +gA) 2 a(p 3 ) and d$gA = -Ag 2 n q ^— / dp 3 

16vr 2 4tt J_ 00 T(p 3 l B ) 

where ^gn q denotes the current of the massless quarks. The initial conditions are given by 



n,(t = 0) = 0, a(ps,t = 0) = l, d a(p 3 ,t = 0) = J ^J^f^ , A(t = 0) = and d A{t = 0) = E . (24) 

y u(p 3 i B ) 

By solving these equations we can see how the electric field vanishes owing to the production of the NO modes. 
Furthermore, we can obtain the temporal behaviors of the electric current densities of the quarks J q = -ign q and the 

gluons J g = f—x '^(pIi'b)^ ' Both of them vanish at t — 0. After the electric field is switched on at t = 0, 

the pair production of the quarks arises and their electric current flows. Similarly, the gluons as the NO modes are 
produced in the vacuum fluctuations and their electric current flows along the electric field. Owing to the production 
of the quarks and the gluons, the electric field decreases and vanishes at t = t c > 0. Hence we compare the electric 
current of the quarks with that of the gluons at the time t = t c when the electric field vanishes, 



gy/gB r+oc , (p3+gA(t=t c )) a(p 3 ,t=t c ) 2 

R(t = t c ) = ^= 4 - y*> , (25) 

Jq Agn q (t = t c ) 

where all of the integration range of the momentum p 3 is taken. But because the relevant modes we should take 
into account are the NO unstable modes, the integration range should be limited to the range in which each mode 
a (P3,t) passes the unstable period until the electric field vanishes. The range is approximately given such that 
\p 3 + gEot c \ < \fg~B for p 3 < and p 3 < \fg~B for p 3 > 0. As our analysis of the gluons is classical, our estimation 
of the ratio is not rigorous. Thus, rigorously speaking, we do not know the appropriate integration range: If we take 
all of the integration range, we take into account irrelevant modes to NO instability, which do not pass the unstable 
period. Therefore, we may take the integration range | j>3 1 < \fg~B in order to see roughly how large amount of the 
gluons are produced compared with that of the quarks. 

In Fig.l we show the temporal behaviors of the electric field and the ratio R with the parameters g = B = Eq = 1. 
We have taken the integration range \p 3 1 < y 'gB and checked that the result does not change even if the integration 
range \p 3 \ < \.2\JgB is taken. We find that when the electric field vanishes, the electric current of the gluons produced 
is approximately eighty times larger than that of the quarks. Owing to this fact, the life time of the electric field is 
much shorter than the life time in the case of the decay only by the production of the quarks. Actually, the above 
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FIG. 1: temporal behaviors of electric field 10E(t) (dashed line) and the ratio R(t) (solid line) 

equations can be explicitly solved when the contribution of the gluons is neglected, i.e. a = 0. The solution of the 
electric field is given by 



E = E cos(^l). (26) 

The solution represents a plasma oscillation 0, HJ. Hence, the life time when E vanishes is given by n 2 /^/g^B, which 
is equal to tt 2 ~ 10 with g = B = 1. It is roughly 8 times longer than the life time shown in the figure. The life times 
are given by 1.3QJ 1 and IOQjT 1 respectively in the physical unit Qj 1 ; Q s ( = lGcV ~ 2GeV ) denotes saturation 
momentum of high energy heavy ion collisions in RHIC or LHC. In this way the decay of the electric field is mainly 
caused by the anomalous gluon production, that is, the production of the NO unstable modes. The contribution of 
the quarks is negligible. 

The NO modes are generated by the vacuum fluctuations and are amplified by the magnetic field. At the same 
time, they are accelerated by the electric field. In Fig. 2 we show the temporal behavior of the momentum distribution 
(j) p or a(ps,t)/ \/T{p^Ib), that is, the growth of the amplitude with time. We can see that the vacuum fluctuations 
( shown by the curve at t = in Fig. 2 ) give the momentum distribution symmetric in p% with the peak at p^ = 0. 
The peak moves to points with negative momentum with time. This is because the gluons are accelerated by the 
electric field so that the momentum with which the amplitude a(p^,t) / \JT{p^Ib) has the largest growth rate is given 
by P3 = - fo dt'gE{t') = -gA(t) < 0. 
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FIG. 2: amplitudes a(p3,t)/ \/T(j>zIb) at t = ( solid line ), t = 0.5 ( dashed line ) and t = 1.0 ( short dashed line) 

We can see, however, from the figure that the peak is not located at p% = — f Q dt'gE(t'), but at momentum near 

P3 = as long as t is small. This is because the initial condition ^"^''"^ = \ TJ , 9B , ^ becomes larger, as I03I 

becomes larger. Thus, as long as t small, the initial condition dominates the growth rate of the modes. Indeed, 
at t < 0.5 the maximum of <p p oc a{pz,t)J \JT(pzIb) is given by 4> p (p3 ~ 0). However, it approaches the value 
4> p (p3 = — Jq dt'gE(t')) as t goes beyond 0.5. The effect of the initial condition diminishes with time. 
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V. GLUON CONDENSATION 

We have mentioned the analogy between the NO unstable modes and unstable modes in the model of the complex 
scalar held. The model describes Cooper pairs in superconductors. Obviously, the unstable modes with zero momen- 
tum in the model form a Bose condensate (ip) = U m 2 /A of the Cooper pairs. That is, once the modes are produced, 
the amplitudes of the modes exponentially grow to the value of the stable state. They are stabilized by the nonlinear 
interaction A|^| 4 . Hence, it is natural to expect from the similarity that the gluons of the NO modes may form a 
Bose condensate [l3j . 

The NO modes typically carry the momentum p = (jot — \J gB,p% ~ ^/gB). This is because the modes have the 
spatially transverse extension given by 1/^/gB and their longitudinal momentum p 3 given by Ap 3 = —gEAt with the 
acceleration by the electric field within the time At = ^/gB/gE. Among them, the mode with the vanishing transverse 
and longitudinal momenta ( pr — and pz + g A = ), has the largest amplitude and growth rate. Furthermore, 
the amplitude of the mode can grow to of the order of 1/g until the nonlinear interactions among the NO modes 
are operative, as we will show below. Hence, we may think that the mode form a Bose condensate of gluons with 
zero momentum. But, contrary to the expectation from the similarity to the model of the complex scalar field, the 
nonlinear interactions in the gauge theory do not stabilize the exponential increase of the mode. They are three 
and four point functions of the gauge fields. The interactions cause the momentum transfer from the mode to the 
other stable modes in higher Landau levels or they produce new NO modes under new magnetic fields produced 
by the electric current of the old NO modes, etc. ( see our comment below ). These new modes induced by the 
interactions have larger momenta than those of the old NO modes. That is, the interactions make the modes with 
higher momenta rapidly increase. In general the number density of the gluons is proportional to the square of the 
amplitude of the gluon field <fi. Hence, once the nonlinear interactions are operative, the number of gluons with larger 
momenta increases, while the number of the gluons in the condensate decreases. In this process the number of gluons 
is not conserved. 

We can estimate how large the mode with zero momentum grows until the nonlinear interactions are operative. 
Namely, we estimate the maximum number of the gluons with zero momentum in the condensate. We note that the 
nonlinear interactions are given schematically by <? 2 $ 4 or gA<& 2 in eq([T]); the background gauge fields A are assumed 
to be of the order of 1/g. They are smaller than the kinetic terms of <I> and A^ when the amplitude <£> is much less 
than of the order of 1/g. But when 4> ~ 1/g, all of the interaction terms become of the same order of the magnitude 
as the kinetic terms. Hence, the amplitude of the mode with zero momentum can grow to of the order of 1/g, in other 
words, the number density of the gluons becomes of the order of 1/g 2 . The most of gluons composing the background 
electric field become the gluons in the condensate. In this way, the gluon condensate with zero momentum arises on 
the way to the thcrmaliztion of QGP. 

Here we wish to make a comment on the numerical simulation 7] by Berges et al.. We can see in the simulation that 
after the NO modes as the primary instability grow large, the secondary instability arises. This secondary instability 
has been discussed[9j to arise owing to the excitation of new NO unstable modes under azimuthal magnetic field, 
which is produced by the electric current J carried by the NO modes. Namely, the electric current J of the NO modes 
becomes large as the modes grow so that the azimuthal magnetic field Bt induced by J becomes strong. When Bt 
becomes strong enough, new NO unstable modes are induced. These new NO modes lead to the secondary instability 
and possess larger momenta ( ps > \fg~B ) than momenta ( pz < V gB ) of the NO modes leading to the primary 
instability. Obviously, the secondary instability arises owing to the nonlinear interactions Q; Thus, the appearance 
of the secondary instability would be onset of the thermalization of the gluons. Actually, the growth rate of the 
secondary instability is much larger than that of the primary instability. Furthermore, the growth rate becomes larger 
as the momentum | ^03 1 of the modes is larger, while the growth rate of the primary instability becomes larger as |p3| is 
smaller. These facts implies that after the gluon condensate with the number density of the gluons of the order of 1/g 
is formed, the rapid momentum transfer arises from the modes with small momentum ( < \fg~B ) to the modes with 
large momentum ( > y 1 gB ). That is, the number of the gluons with the large momenta rapidly increases owing to 
the effects of the nonlinear interactions. The process would lead to the equipartition of the momentum and eventually 
thcmalization of QGP. This simulation indicates that the gluon condensate may arise before the thermalization. 

It is interesting to speculate that inverse cascades, with which the modes with small momenta are amplified in a 
model of scalar fields shown in the recent paper[I3], correspond to the excitations of the NO modes and cascades, 
with which the modes with large momenta are amplified, correspond to the excitation of the modes of the secondary 
instability. We need to analyze more closely the correspondence in order to confirm this speculation. 
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VI. SUMMARY AND DISCUSSION 

Motivated by the recent study of the anomalous production of the Nielsen- Olesen unstable modes by the Schwinger 
mechanism, we have discussed the decay of the color electric field in the classical approximation by taking into account 
the back reaction of the electric field to the gluon and quark production. We have found that the electric field rapidly 
decays owing to the anomalous production of the gluons. The contribution of the quarks to the decay is negligible. 
We have also found that the amount of the produced gluons is about a hundred times larger than that of the quarks. 

A model of color glass condensate predicts that the color electric and magnetic fields are produced immediately 
after high energy heavy ion collisions. Fluid dynamical simulations of thermalized QGP suggest that they should 
decay into the plasma within a time 1 fm/c. Our analysis indicates that such a very fast decay is caused by the 
anomalous gluon production. Actually, our analysis shows that the decay is completed within the time of the order 
ofQ-\ 

We have also discussed a possibility that the gluons of the NO unstable mode form a Bose condensate with zero 
momentum until the nonlinear interactions in the gauge theory are operative. The number of the gluons in the 
condensate rapidly increases and then saturated when it becomes of the order of 1/g 2 . After the condensate with such 
a large number of the gluons is formed, the nonlinear interactions lead to the rapid increase of the number of gluons 
with large momentum ( > y ' gB ). Hence the equipartition of the momentum and eventually the thermalization of 
QGP would be achieved. 

We have discussed the decay of the color electric field in glasma. The color magnetic field in glasma also decays 
in the following. In general the longitudinal color magnetic fields form flux tubes, which expands with time. Owing 
to the expansion, electric field SEt perpendicular to the magnetic field is induced according to the Faraday's law of 
induction. On the other hand, owing to the expansion of the longitudinal electric flux tubes, a magnetic field SBt 
is induced, which is parallel or anti-parallel to the electric field 5Et- Thus, under the field 6Bt, NO modes 5<f>T are 
excited and make the electric field SEt decay rapidly. Eventually, the expansion of the magnetic flux tube induces 
the electric field 5Et, which decays by the acceleration of the NO modes 84>t- This implies the decay of the magnetic 
flux tube. This is the decay mechanism of the magnetic field. 

We have discussed the gluon production in the non-expanding glasma. When we treat it in the expanding glasma, 
the similar analysis is possible. But, we should use the initial conditions for the NO modes shown in the recent 
paper 15]. Using the initial conditions, we can properly take into account quantum effects of the modes under the 
background gauge fields. 



We would like to express thanks to Dr. K. Fukushima for making me pay attention to the recent paper by Berges 
and Sexty. 
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